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SUMMARY 


An unsteady lifting surface theory for the counter-rotating propeller Is 
presented using the linearized governing equations for the acceleration poten- 
tial and representing the blades by a surface distribution of pulsating acous- 
tic dipoles distributed according to a modified Blrnbaum series. The Blrnbaum 
series coefficients are determined by satisfying the surface tangency boundary 
conditions on the front and rear propeller- blades . Expressions for the com- 
bined acoustic resonance modes of the front prop, the rear prop and the combi- 
nation are also given. 
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INTRODUCTION 


Interest In propeller propulsion has been revived recently by the need to 
' economize on fuel consumption and fuel costs. Besides, with the development 
of transonic airfoil sections, even with airscrew propulsion. It Is now pos- 
sible to attain cruise speeds comparable to those of current turbofan trans- 
port aircraft at similar operating altitudes. Preliminary studies appear to 
Indicate possibilities of realizing a clear 15 to 20 percent superiority (see 
Mlkkelson et al.) of propulsive efficiency for the advanced turboprops over 
that of turbofan propulsion. For commercial acceptability It Is Important to 
keep the noise levels of the propeller system low. Consequently, NASA has 

I 

mounted a program for advanced turbopropeller development to realize these 
objectives . 

It has been known for a long time that both the power absorption charac- 
teristics and the efficiency of propellers at high forward speeds may be 
Improved by Increasing their solidity. Solidity can be Increased by Increas- 
ing the blade chord and the number of blades. In addition, the use of a 
counter-rotating propeller mounted closely behind on a common axis provides 
great Increases In the efficiency since the energy losses due to slipstream 
rotation can be eliminated. When the two propellers are close together both 
of them will have nearly equal Inflow velocity. But 1 , due to the rotational 
component Imparted by the front propeller, the rear propeller operates at an 
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effectively higher rotational, velocity or equivalently lower advanced ratio 
than the front one. 

An extensive survey of recent NASA theoretical and experimental research 
on propellers has been published by Mlkkelson (ref. 2), Bober and Mitchell 
(ref. 3). Methods available for the aerodynamic analysis of propellers may be 
classified broadly Into Euler equation methods (ref. 4), vortex lifting line 
theory (refs. 5 to 11), lifting surface methods (refs. 11 to 18) and the Ffowcs 
W1 1 1 lams-Hawklngs equation solution method (ref. 19). It will be seen that a 
calculation of the acoustic Intensity distribution In the flow field of the 
propeller requires a knowledge of the aerodynamic loading on the propeller 
blades. Hence It Is necessary to determine the pressure distribution on the 
blade surfaces prior to a noise field calculation. 

A lifting surface theory will be presented here for a subsonic counter- 
rotating propeller (fig. 1(b)) using the acceleration potential formulation. 
Unlike the velocity potential perturbation, this method has the advantage that 
It does not have to deal with the discontinuities arising from the potential 
vortex sheets emanating from the front and rear propeller and their Inter- 
action with the blades. It will enable us to determine both the aerodynamic 
load distribution on the blade surface and the acoustic pressures In the flow 
field concurrently. We shall adopt a coordinate system In which the front 
propeller Is stationary so that the free stream approaches It with an axial 
velocity W^ and a rotational velocity n^r. In this coordinate system, the 
rear propeller rotates with an angular velocity, ± (+ sign for counter- 

rotation and - sign for corotation). .The axial velocity Is assumed to be the 

' / . ' 

same as for the front propeller but the airflow acquires a circumferential 
velocity (Q^r - W fl tan e^) In passing through the front propeller. The rear 
propeller removes this rotation and converts it Into useful thrust. The sur- 
face boundary conditions on the front propeller blade are satisfied In a 
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rotating coordinate system while on the rearl propeller they are satisfied In 
the same moving coordinate system but with the blade surfaces nonstationary. 

This method gives a unified treatment of both the front and the rear pro- 
peller of a dual rotation propeller system taking Into account their mutual 
Interference In unsteady subsonic axial flow. Thus, It Is possible to calcu- 
late the resultant aerodynamic loading on the blade, the performance of the 
dual rotation propeller and Its acoustic field at any point In the flow for 
blades of arbitrary geometry. Expressions will be given for the blade surface 
loading, the acoustic pressure at a point In the flow and the acoustic reso- 
nance characteristics for the front propeller, the rear propeller and the 
combination. 

FORMULATION 

We now consider a counter-rotating propeller with the angular velocities 
and for the leading and the trailing propellers of equal hub and tip 

diameters, W the axial velocity of the airflow Into the leading propeller. 

The leading propeller develops a thrust T^ by accelerating the free stream 

by the Induced velocity w^ In the wake far downstream. In addition, the 
propeller also Imparts a rotational velocity to the airflow In the same direc- 
tion as . When the number of blades In the propeller Is large, the 

air passing between the blades Is guided closely so that we may assume the 
flow through It Is similar to that through an axial compressor cascade except 
for the different boundary conditions at the tip. The air angle c* 2 at the 
trailing edge nearly equals that of the mean line of the passage between the 
two blades. Assuming that the two propellers are close together, the axial 

i 

velocity for both of them may be regarded as nearly equal. The free stream 
velocity components for the front and the rear propellers are given In cylin- 
drical coordinates by 
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(0,fi.|r,W a ) and (O.fi-jr - W fl tan <* 2 ,W a^ 

In this event, the linearized differential equations for the pressure field of 
flow singularities on the blades of the counter rotating propeller are the same 
as those for the blade rows of a stage of an axial compressor. 

We shall represent the propeller blade surface by pressure singularities. 
Since the rear propeller lies downstream of the front propeller, It experiences 
a periodic wake which would cause the pressure pole on the rear propeller blade 
to fluctuate In strength with time. When the free stream has periodic fluctu- 
ations, of angular frequencies and u 2 for the front and rear propellers, 
the nondlmenslonal linearized differential equations for the perturbation field 
of an oscillating unit pressure monopole on a blade of the front and rear pro- 


peller situated at p-j (p^, , C-.) may be written respectively as 


3 P 1 1 3p 1 /I 

2 r, ar, 2 

ar i 1 1 Vi 
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In which we have defined the dimensionless parameters 


r l 

= r/R 

Z-, - z/R 

X 1 

= W /Q. R 

3 1 

H 2 = M tan o 2 

P 1 

= p/R 

?! = ?/R 

X 2 

= W a /n 2 R 

B =Vl - H 2 

l l 

= tn o 

t = R/W 
o a 

S2 = 


B 2 "V 1 ' M 2 

W 1 

= OJ-j /fi-| 

u = w 2 /q^ 

fi r 

= fi-j + « 2 



( 2 . 2 ) 

i > 

Since the rear propeller Is rotating with an angular velocity (« 1 + Q 2 ) rela- 
tive to the front propeller, the azimuth angle cp of the pressure pole Is also 
changing with time as (q> + fi r t). We express the perturbation pressures p^ 
and p^ In terms of the Green's functions Pi ( r ^ and P^r^, 
p-l ; k 2 , a 2 ) as 


J fCO 

l[k 

P l ( »P"| »^i ® 

0 

J r* co 

P 2( »P"j * ^2 » ^ 

o 


( 6-<p) +a^ ( z-j -C-j ) t^ /X-j ] 


da n 


(2.3) 


1[k 2 (e-q>)+a 2 (z 1 -c 1 )+o>t- ) /\ 1 ] 


da. 


(2.4) 

Combining equations (2.1), (2.3), and (2.4) and equating the corresponding 
terms, we get the differential equation for the front and' rear propeller mono- 


pole Green's functions P-j and given by 


i_ ffi 

dr* * r l dr l, * 


21 1 * "if 2 k ? 


"'VI 


- a 


1 


1 2 
r l 


P 1 - " 


4(r 1 - p 


dr ? V' 4 r 2 r? ! ‘ ' AJ' ! 


(2-5) 
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where 



P 2 = P 2 r ] /2 * “ = 0 2 k 2 

u + k 2 0 + S) 

“k = h 


MM 2 k 2 (a 2 + w k ) 



( 2 . 6 ) 


We can expand the delta functions 6(6 - q> + wt^/X^) and A(z^ - f^) In 
Fourier series-integral form as 


6(e - «p)6(z. 


(«-v) + 1a 1 ( z, -C-, ) „ 


6(0 - tp + wt-j /X-j ) 6 ( Z, - f i ) = 2 ^ / e 

4ir k 2 */-“> 


1k 2 (0-q>)+1a 2 (z 1 -C 1 ) 


(2.7) 


which can be used to determine the Green's functions P^ and P^ for the 
monopoles on the front and rear propellers governed by the nonhomogeneous 


differential equations (2.5). 
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FLOW SINGULARITIES FOR THE TWO PROPELLERS 
Equation (2.5a) for the front propeller Is a nonhomogeneous Bessel's 
differential equation while equation (2.5b) Is the nonhomogeneous Whittaker’s 
differential equation. We can write the complementary solutions for these two 
equations as 

VVV = A*J |< ^(X k ^r 1 ) + B * Y k . | ( x k . i r -|) 


wv ■ ) c *\ 2 . s ( Vi ) 4 M v^ , V' ) j’ k V' /2 


(3.1) 


These solutions must satisfy the mixed type boundary conditions given by 
P-, = 0 at r ] = 1 p 2 = 0 at r -| = R 12 

dP dP 

" ° at r l - h l " 0 at r l ■ R 12 h 2 

(3.2) 


where A*, B*. C*, and D* are arbitrary constants. We shall assume here 
R-J 2 = ratio of tip diameters of the front and rear props = 1. We look for 
solutions such that there Is no radial flow at the hub and the pressure per- 
turbation vanishes at the blade tip. From these boundary conditions we can 
obtain a transcendental equation 


J k 1 -l (x k 1 h l ) " J k 1 +l (X k 1 h l ) J k 1 (X k 1 ) 

Y k 2 -l (x k 1 h l ) ' Y k 1 +l (x k 1 h l ) 

giving the eigenvalues x, ,8. for the front propeller. The corresponding 

K 1 

Green's function for the oscillating unit pressure monopole on the front 
propeller Is given by 


(3.3) 


Y k ^ X k ^ 

= ' \(\ ) J k 1 (x k 1 r l ) + Y k 1 (x k 1 r l ) 


(3.4) 


8 


The Erdelyl's form of Whittaker function used In equation (3.1) can be 
related to generalized Laguerre functions through Kummer's function by the 
relations 


N 


z a-l/2 

X,a (Z) = r(l + 2a) 


M 

X.« 


(z) 


e- z/2 z 2a . . 1 

r(l + 2a) l F l (a + 2 ' x! 


+ 2«; z) 


-Z/2 2a /x + a 



(3.5) 


where M Is Whittaker's first function, ,F. Is Kummer's confluent hyper- 

x,a 1 1 

geometric function. In general there are three distance cases giving 

transcendental equations for the eigenvalues of the rear propeller. We shall 

consider here only the case (1 ± 2a) * - m, m = 0, 1 , 2, 3, ... for which 

the eigenvalues are given by 


5 , <2») . (2i) 

b l - 2 V b loU, < h 2 y) - b 3 b 4 L— < h 2 y > 

J_V l b B 

5 . (-2S) , (-2S) 

b 2 " 4 h 2 y b 3^— ? * h 2 y * " b 4 b 6L <h 2 y) 

5 b 3 I b 7 

In which 


. (2a) 

b ioL (y > 

5 b 10 


5 b 3 



(3.6) 


(3.7) 
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The Green's function for the pressure monopole on the rear propeller Is there- 
fore given by 


P 


2 




1 r 
2 2 



where the coefficients (C,D) are given by 


(3.8) 


L \ «-**/ 

(n) 

F (f b l?) f b 4 

r (f b n) (n) 

2 . 

5 b 10 



(3.9) 


It Is seen that the pressure monopole on the propeller has a characteristic 
Bessel function variation radially. The pressure monopole on the rear propel- 
ler has a characteristic radial variation described by the Whittaker's func- 
tion. In the presence of an external forcing function, the differential 
equations (2.5) are nonhomogeneous and the pressure variation due to the mono- 
pole may be expressed by expanding In terms of appropriate orthonormal eigen- 
functions. We shall choose the eigenfunction for the front propeller as the 
function $ (r ) given by a linear combination of Bessel functions of the 

1 *» I 

first and second king, J n (r-|) and Y n ( r -|) normalized to satisfy the orthogo- 
nality condition 



r l P I ( Nci r l )P l (x km r l )dr l * 


<$_ m. „ 
SLm kil 


w 


P l (X ksi r l )/m ksi 


(3.10) 


For the second propeller, since the Whittaker functions do not possess orthog- 
onality properties. It Is convenient to relate them to generalized Laguerre 
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polynomials using the relations of equation (3.5). The normalized generalized 
Laguerre polynomial 1 ^ 2< ^(r 9 ) satisfies the orthogonality condition 


12 . (2 ‘> 


R 1 2 h 2 


|__ (r 2 ) |_ (r 2 ) dr 2 = n 


am aka 


, (2a) / 

( r -. ) = I — _ (r ? ) / n_ 

aka 8. c / ak 


The solutions of the nonhomogeneous equations (2.5) may be expressed as the 
eigenfunction expansions given by 

P 1 (W a i) = l 9 a 1 5 ( a i )^k n 9. ( r l )# k ] a (p l J 

(k„) (k ) a+ } 

P (r P ;a ) = 2 £ h (a )*_ (r )^ (p,)(r P ) 

2 1 1 r * 11 2 lk4 1 22 

- I (r 2* p 2 ) r: 1/2 » ,/2 
e 2 ik », 


<k.) (k I 

where g (a ) and h (a 9 ) are coefficients of the expansion to be 

U 1 

determined. We shall now expand the Dirac delta functions 6 ( r^ - p^) In 
equations (2.5) In terms of the eigen functions corresponding to the front and 


rear propellers as 


r i V " yVi *k 1 a (p i ) *k 1 a (r i ) 


~ p ( p + ^p ) z 

Hr, -P,) =He (r,P,) a *_ ( P ) 

1 1 ~ a £ £ ak 2 a 1 


(3.13) 


n 



Then, the eigenfunctions §. (r.) and ¥ 

K- r 1 ak 2 a 

the Bessel and Laguerre differential equations 


(ri) 


satisfy respectively 


d 2 # 


dr 



d$ 


1_ 

r. 


k-j SL 


dr. 




2 

d? 


8,k 2 8. 


dri 


(V-' ) 


dS^ 

aTkpB. ^ 

-j — - + 7 - f 

dr 2 r 2 ak 2 a 


(3.14) 


where a Is the degree of the Laguerre function defined by the relation 

-3/2 


I - x 


ak 2 a 


ak 2 a 

Tk^7 

aa 


-(=♦« 


(3.15) 


in terms of Its eigenvalues v . The Independent variables p and r 

ak 2 a 1 l 

are defined In terms of r^ by the relation given In equation (2.6). In the 

above are the radial eigenvalues of the front and rear propeller, 

respectively. The elgensolutlons of the Laguerre differential equation exist 

If the degree a Is an Integer. Considering the oscillating unit pressure 

(k ) 

pole on the front propeller, the expansion coefficient g^ Is obtained by 
substituting for and 6 ( r ^ - p^ ) from equations (3.12a) and (3.13a) Into 
equation (2.5) and equating the corresponding terms. Thus, 


where 



(a ] ) 





M = M(1 + k ] )/X 1 


(3.16) 


(3.17) 
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Both M and x k are always real. When the axial flow Is subsonic, M < 1, 
2 ^2 

0 > 0 and the sign of A k ^ depends on the relative magnitudes of X k 

and M/0. But, for supersonic axial flow M > 1, B 2 < 0 and A k Is 
always positive. 

The perturbation pressure p^ In equation (2.3a) may be expressed In 
terms of the eigenfunction $ k a using the Fourler-Bessel expansion (3.12a) 

as 


_ . , dO - 1[k (e-q>)*» t /x_ ] 

P i< r i* p rt> ' 72 l <V e 

4ir 9. 1 1 


(3.18) 


where 


(k 


) „ r la Z (k ) 

(Z-,) = I e g^ (a 1 )da 1 , Z-, = z., - C-, 

V O 


(3.19) 


The Integral may be evaluated by contour Integration methods and written as 


(k ? ) - 

X 9 (Z 1 )= X^0 ' e 


Z 1 (l MM/02 - A k j/ej Z^l MM/0 o +A l , ft /0 


2 k.j 9. 


(3.20) 


If A k ^ > 0, the second term In equation (3.20) Is rejected since the 
Integral diverges with Z^ . We therefore put 


/Vj , _ /Vi 

* V 




(3.21) 
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Thus, the solution 

Is always bounded 

and the sign of 


Is assigned as 

follows : 





If 

A k lt > 0 

put = 

•V 


If 

A £ a < 0 . M > 0 

= 


1 

If 

< 0 , M < ° 

= 

’'V 

1 sgn Z 1 


(3.22) 

Thus, the pressure field of a pulsating unit pressure pole situated at 
(p, , <p, f-j) on a propeller blade becomes 


M k, l ** 


r i > *k 1 a. ( p i 5 A k n a. e 


1[k 1 (©-(p)+f 0 Z 1 + Ul t 1 /\ } ] 


The field of a pressure pole situated at (p^ , cp + <p + , C-|) on the 

-th blade of the front prop Is given by 


= 4irB l j[ *k 1 l <r l ) *k 1 l tp l )A k 1 l 


1 [ k i ( 6-<P-| « 1 4 1 +f 0 Z 1 t-j /X-| ] 


(3.24) 


where Is the uniform Interblade phase angle factor of the front prop, 
a . j = 0, 1, 2, 3, ..., (Zi - 1), and 6^ Is the Interblade azimuth pitch angle 
2ir/z.| . The pressure field of a pulsating unit pressure dipole situated at 
p-j , ip + <p^ + m.^, ^ ) Is obtained from that of the pole by differentiating p^ 
along the normal to the local camber line giving 
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(m ) 1 cos a, , 

*D = " 4ir0 £ J ^k-j r i a (p l ^ A k n a 


x f 2 e 


1 [ k 1 ( Q-^-, -<p) +m n 4-j ( o 1 ) +f Q Z 1 + Wl t 1 / v, ] 


Putting 


1 

f 0 = — - f tan a, 
2 p-, o 1 


(3.25) 


k l = Vl + 5 1 


(3.26) 


the pressure field of unit dipoles situated on all the Z^ , blades of the front 
propeller Is obtained by summing over m n and we get 

"L > 


Z l _1 (m, ) 1Z, COS a, 

P |3 = E Pj) = ~ 4^3 ^ ^ *k 8^ r l^*k st^ p l^ A k a. 

U 1 m=l u 4 * 13 k 1 a V 1 V ' V 


x f 2 e 


1 C k i ( 0-q>-| -<p) +f qZ^ +u>-| t-j /\-| ] 


(3.27) 


Again considering an oscillating unit pressure pole on the rear propeller, the 

(k ) 

expansion coefficient h c (a _) Is obtained by substituting for P ? and 

U £ 

5 ( r^ - p^) from equations (3.12b) and (3.13b) Into equation (2.5) and equating 
the corresponding terms as 
(k > 

<a 2 > = 

n 


r 2 

«fi 3 [a r 2HM 2 k 2 (a 2 + £ k )] [a 2 - a 2 3 2 ) - a 2 


3/2 


M 2 ' 2 

M o, k 

0 = 3 


4 • l 1 = l + o + 2 


(3.28) 
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The perturbation pressure p 2 In equation (2.3b) may be expressed In terms 

of the eigenfunction ¥ (rj using the Fourler-Laguerre expansion 

aKpIl 

equation (3.12b) as c 


p2^ *"i * P-| * ) 


1 1 . . V 

1 v v v 1/2 / \ iti , w a 2 ' 2 ^ r 2 p 2' 

— j £ I I V (r,)*^ (p n )r ? P? e 
4v ~ k 2 t ak 2 £ £k 2 a fck 2 & 

(k ? ) . 

(Z 2 )exp[1k 2 (e - <p) + 1(c + kgQ)^/^] 


(3.29) 


where 


(kp) . r la ? Zp (k ) 

(Z 2 ) = I e * X, * (a 2 ) da 2 Z 2 = z - f (3.30) 

ii / n 

«/ —co 


This Integral may be evaluated by contour Integration methods and written as 


where 


(k„) . , 1f,Z„ 

Jj (Z 2 ) =-^e 42 S, 

ll 1 68 M 


S* = 


3/2 fl o = 2MM 2 k 2 0 2 


, -2 

L 1 “k 


H 2 k 2 ( 4 o ' a o) 

f *' !? 0 * *4 S9 " z ) 


(3.31) 


(3.32) 


We can then write p 2 as 

- 1 1 / % 

. . * 1/9 Z a '*’ 9 “ 9 V^o^Pp) 

P 2 (r ,p ,t, ) = — L -j- l l <r )*. (p,)r“p d e l 

32*8 M ~ k 2 a. ik 2 a ak 2 s. ak 2 i 


1f 4 Z 2 +lk 2 (e-«p)+1(w+k 2 fi)t 1 /\ 1 
S*e (3.33) 


giving the field of a pulsating unit pressure pole situated at the point 
(p^, <p. C 2 ) on the rear propeller. When the pole Is located at the point 
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(p-|, <p 2 + q> 2 + m 2 6 2’ w1th a un ^ orm Interblade phase angle o 2> the 
pressure field Is given by 



”3 l * 
32ir0 n ~ k ? 


I V I /2 ('i)*. (P n )r?P ? 
a ak 2 a j.k 2 i s.k 2 a 


(r 2 +p 2 ) 


S*exp 1 f 4 Z 2 + k 2 (e - <p 2 - <p - m 2 6 2 ) + m 2 o 2 6 2 + (u + k 2 Q)t 1 /x 1 

(3.34) 


Again, the pressure field of a pulsating unit pressure dipole at the same 
point on the blade Is obtained by differentiating p 2 normal to the blade 
mean camber line and Is given by 


( m 2 >-. * . . ~ C0S “2 . ■■ „ 1/2 


2 I »_ *_ <r,>*_ (p,)r * p, 

I 1 1 u .. b « a ft 1 ft a 1 c c 


- 1 1 / 

<*+ ~ - « (r 2 +p 2 ) 


2 e 2 


32ir|3 M ~ k 2 8. «.k 2 9. fck 2 & 


a.k 2 ». 


[< 


f <; 6Xp + “ V9 _ < p ) ~ m 9 6 9 ( k 9 - a 9 ) + 


2 ^ 2 ' 2 


(w + k^Jt-j/VjJ 


where 


- — + f. tan a 0 )S* 
p 1 4 * 


(3.35) 


(3.36) 


& 2 = 2tr/Z 2 = Interblade azimuth pitch angle of the rear propeller and = 

Interblade phase angle factor of the rear propeller, = 0, 1, 2, 3 

(Z 2 - 1). For the Z 2 blades, the perturbation pressure of the unit dipole 
situated on all the blades of the rear propeller Is obtained by summation over 


m 2 and Is given by 
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P D 2 (r l ,P 1 *V 


V 1 (m ? ) _ , 

J , Pd, ‘V'i-V 

m 2 =l 2 


Z 2 cos « 2 


1/2 


rH I *. (r,)*, (pjr“ p 

I. A Al. n Al, n 1 ft L, ft 1 c 


- 1 1 . , 
- a + ^ - 2 < r 2 +p 2> 


32ir6 M ~ k 2 a fck 2 & ak 2 t a.k 2 a. 
f 5 exp l[f 4 Z 2 + k 2 < e - <P2 " ^ + + k 2 




(3.37) 


where 


k 


2 


U 2 Z 2 + °2 


(3.38) 


The dipole solution for the second prop consists of an axial component which 

Is wavellke In Z, the axial distance from the rear propeller. 

2 2 

When (A q - M*) -» 0, S* -* <» and p Q -» ® and Is termed acoustic 
resonance. At resonance, A q = ±M*. The eigen-numbers a and k 2 at 
resonance satisfy the relation 


8. = 


2 ~ ®2 k 2 + 


2H 2 M 2 k 2 

B 2 




(3.39) 


the + and - sign corresponding to the advancing and retreating waves 
respectively. 

In equation (3.27) If A. = 0, p n -> ® and we have acoustic resonance 

In i SC L/-| 

of the front propeller. The corresponding circumferential mode number and the 
eigenvalue x k ^ satisfy the relation 


M(1 + k^ 
Vjl “ 4 Bx^ 


(3.40) 


The resonance condition for the counter-rotating propeller system Is obtained 
by setting the azimuthal eigen-numbers k^ and k 2 of the front and rear 
propeller equal. We then obtain a relation between ? and x^ ^ given by 
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I = - \ - 0 2 (-l ± — 


(-. . . .. 


(3.41) 


The condition k leads to the relation 


°1 = « 1 Z 1 - v 2 Z 2 


(3.42) 


connecting the Interblade phase angles of the front and rear props. 

Equations (3.27) and (3.37) give the pressure field of a pulsating unit 
pressure dipole situated at a point of the front and rear propeller blade. 

If Ap.j and Ap 2 be the upward acting pressures at the point (p , <p, C-|) and 
(p-j. <p. C 2 ) of the front and rear prop blades due to both the camber and 
thickness distribution and ^ , a 2 are the blade angles of the front and rear 
propellers at radius p^ , considering an element of area (dp^ dc^/cos ) and 
(dp^ d£ /cos a 2 ) at these points, the strength of the pressure dipole at these 
points would be 

front: Ap^ dp^ d^/cos rear: Ap ? dp^ df^/cos « 2 (3.43) 


The perturbation pressure field P-j and P 2 for the whole blade Is given by 

1Z 1 f 1 f ZlTE 

1 4«B J J ap 1 (p 1 ,<p,f 1 )**(r 1 ,P ^ , ti ) dp^ df^ 


h l Z 1LE 


\ = — V 1 1 i 


32irB M ~ k- !. qjf p . j 

1 2 R 12 h 2 Z 1LE 


/’“ /' 


2TE 


-» 4 


ap 2 (p-| ,<pf£ 2 )'M r i ,p i »*! > dp i d ^i 


(3.44) 


where ( Z 1L £. Z-| T £) and ( Z 2LE* ^2TE^ deno ^ e ^e ax ^ a ^ positions of the leading 
and trailing edge of the front and rear propeller blade at any radius p^ and 
and are the kernel functions defined by 
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-i 1 [ k-i ( ® -t P-i _< p) + f n ^i t, /X-. ] 

*.<VVV . * k] 1 ( r 1 )« kil (P 1 )*i ,f 2 e ollll 




v *< .P-, •t 1 ) -» ( r i>* ( P-j ) r 2 P 


a “ + 2 „ 2 (r 2 +p 2 ) 1/2 


8k ? 8 ik 2 i 


SLk^a 


f 5 e 


1[f 4 Z 2 +k 2 (e-v 2 -<<>) + (u+k 2 n)t 1 /\ 1 ] 


(3.45) 


If the blade loadings Ap^ and Ap 2 are known, the resultant acoustic 
pressure at an arbitrary point In the flow at time t^ Is given by 


P(r^,t^) = P 1 (r 1 ,t 1 ) + P 2 (r 1 ,t 1 ) (3.46) 

Since Is a function of the mode numbers (k, 8) and P 2 Is a function of 

the mode numbers (8., k, 8), the resultant acoustic pressure P at a point Is 
a function of the mode numbers (8, k, 8). Summation over 8 (8 = 1, 2, 3 ...), 

k (k = 0, 1, 2, 3 ...) and 8 (8 = 0, 1, 2, 3 ) of equation (3.46) gives 

the total sound pressure level at any point for all modes together. 
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PERTURBATION VELOCITIES 


We shall split the surface pressure distribution Ap 1 and Ap ? on the 
blades Into thickness and camber components for blade. The thickness and 
camber contributions will be separated Into radial and chordwlse variables as 


4p n - WW 

AP lc = F 2<-l>*|cft<^l> 

np 2t - s,(S 2 )*_^(p,) 

AP 2c = G 2 ( “2 ) \^( p l ) 


where u 


Is a Glauert angle defined by 
y-j = - c i cos <o.| - ^ < y' < + c, 

y 2 = - c 2 cos Z 2 - c 2 < y^ < + c 2 


O < g>i , g > 2 < ir 


(4.1) 


(4.2) 

The functions , F 2 , G 1 , G ? will be expressed In a Blrnbaum-Glauert series 
expansion as 


F-| ( <»>) = cot j +■ l A m sin mu F,(«) = B, 

m=l 


G-j(w) = C Q cot - + l C m sin mw G 2 (u) = D q 

m=l 


(cot r - 2 sin u) + l 

' / m=l 

(cot ^ - 2 sin u) + l 

' / m=l 


B_ sin mu 
m 


D„ sin mu 
m 


(4.3) 

assuming the coefficients A. B. C , D m to be Independent of I, k, and it. 

m m m m r 

Combining equations (4.1) to (4.3), we have 


^1 ( W 1 * ^*1 * ^ ^ ) = 4P "| 
H 2 (u 2 ,p 1 ;l,k,l) = A P 2 


cot f 1 + J^m s1n ^iW^l* 


k cot Y + * s1n (Pi) 

V m=l m 7 ikl 


(4.4) 
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The coefficients and are determined by satisfying the surface tan- 

gency boundary conditions on the blades. To satisfy the boundary conditions 
of flow tangency on the blades, we have to obtain the perturbation velocities 
u and I from the perturbation pressures and P 2 by Integration of 
the unsteady Euler equation of motion. For this purpose we Introduce the heli- 
cal coordinate system (r^, , t^) and (r^, a^, i^) for the front and rear 

prop based on the undisturbed free stream for each (fig. 2(c)). We define the 

helical angles e, and e. defined by 
h l h 2 


tan cL = tan e. - tan 

2 h^ 2 


In addition, we also Introduce the blade section coordinates (r^, y.j , z.j) and 
( r l , y 2 , zp (fig. 28) which are related to the cylindrical coordinates 
(r^, e, z^ ) by the transformation 



(4.6) 


We also have the transformations between the (u lr , u y , up, (u 2r , u^, u^ z ) to 
the (u rl , u^ , u^), and (u r2 , u o2> u^) coordinates 
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(4.8) 


The first matrix transforms from the (r^, , t^) to (r^, y.J , z.j) coordinates; 

the second matrix transforms from ( r^ , a 2 , t 2 ) to (r.j f y 2 , z£); the third 
matrix transforms from the (r^, a 2 , t 2 ) to the (r^, y.j , z.j); and the fourth 
matrix transforms from (r^ , t-j) to (r-j, y£, z£) coordinates. The resultant 

velocity Including the respective free stream components are given by 
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0 

cos e 
sin 0 

0 0 

cos e 2 sin e 

sin e 2 cos e 



( 4 . 9 ) 


The dimensionless equations of motion In the corresponding helical coordinates 
(r, a, t) may be written in the form 

-» -> 

dU-| _ 3U^ 1 fdP 1 3P ] 3P 1 

8*7 * 1,01 ^ ' " y ^ ^ 


( 4 . 10 ) 


3U 2 

3u 2 

1 

f3P 2 

3P 2 

ap 2 

at, 

+ U 02 3o 2 " 

y m 2 

i 3r i 

C\J 

0 

CO 

' 3t 2 
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where 


'01 


( w ! * n i r i ) 

W_ 


1/2 


vl/2 


1 


■ n - sec \ 


02 


W l _ (rft r - W a tan a 2 ) 2 
_ 


1/2 


1 + (^- tan ~ a 2 ) 


1/2 


= sec e t 


so that we have 


u. = u,/W, u- = u./W, 
1 la 2 2 a 


(4.11) 



sec e. 
h i 

-* 

a^ 

l l 

/aP 1 

B A 

3P 1 


3o-| 

” yH 2 

K^\ ’ 

a o 1 ’ 



sec e. 
h 2 

-» 

au 2 

i | 

f aP 2 


8P 2 

V 

aa 2 

2 ' 

V» r i • 

3o 2 * 

3t 2 


(4.12) 


We shall assume that the perturbation velocities u-j and u 2 corresponding 
pressures P 1 and vary harmonically as follows: 


^ 2 Ic^t-j/X-j ^ 2 1(w+k 2 «)t 1 /x 1 

" * u 2 = u 2 e 

1(w+k 2 fl)t 1 /x 1 


U 1 " U 1 e 


P 1 = P 1 e 


1«^ t^ /x^ 


p 2 = p 2 e 


(4.13) 


Substituting Into equation (4.12) we obtain the equations 


au, ^ ■* cos e h. /ap, aP, ap ' 

^ 4 ’ N C0$ e - 1 =1 ■--^“^■^■ 5 ^. 


2 (<*> + kjQ) 


au . 

3o 2 + 1 X 1 


2 cos e h ? /aP 2 aP 2 aP 2 

cos e. u, = - I— , —■ , - — 

h 2 2 yM 2 \ 3r l 3o 2 3 t 2 / 
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(4.14) 



which can be Integrated with respect to and a 2 respectively giving the 
perturbation velocities u^ and u 2 produced by the front and rear propellers 


each In Isolation as 


- S°1 cos 6 h, A i - 
U 1 c ■ V 


tos 6 h /•'“ 1°! C0! 9 h/ x l (it 3P ap 

J e ~ ~ 


ar 1 ’ • 8x^1 


l(u+k 2 fi)/x 1 c» 1 cos e h ^ cos e^. 


2 “ " ,M 2 


■i(«+k 2 «)/x 1 c 2 cos e h a * a p 3 p \ 


(4.15) 

where u^ m and u 2m are the free stream values of the perturbation velocities 
corresponding to ^ , o 2 ■» « which we shall assume to be zero so that we have 


- 1 w_ a, cos e. /X, 

11 h, 1 

e cos e. /'1« 1 o 1 cos e. /X, /. n x 

-» h, / 1 1 h, 1 /3P, 3P 3P, \ 

^ v r 


-1(a>+k 2 Q)/X^ )<? 2 COS e h 


COS 0 L 


•1((0+k 2 fl)/X 1 a 2 COS 0. /gp gp a p \ 


(4.16) 


The helix angles 0. and 0. can be related to the variables 0 and z, 

h.j h 2 1 

by the relations 


z, z,S2 R Qz. 

0. = 0 - r 1 0. =0 - -rj = 0 - — 

h 1 x^ h 2 w a x 1 


(4.17) 
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the differential transformation from (r.j, e, z 1 ) to the helical coordinates 


( r ] , d r t-j) and (r^ a t 2 ) Is given by 


/dr, \ 


r-j de 




h 


V 


l cos e k 
h i 



M 


r 1 de 


v d v 


h 


v° 


V dr A 


0 sin e h -cos 


0 cos e. sin e. 
"2 


da„ 


Jr 


V 


(4.18) 


When moving on a given helical stream line, = constant, d^ = - 0 and 

we have 


front prop: dc^ = dz^ sec ; rear prop: do 2 = dz 1 sec 

Integrating and putting z = 0 for a = 0, we have 


(4.19) 


cr. = z, sec e. 

1 1 h. 


cr_ = z. sec e. 

2 1 h. 


(4.20) 


Hence, performing the Integrations In equation (4.16) along the helical 

streamlines of the respective free streams for the perturbation velocity 

i 

vectors we obtain 



(4.21) 

where the amplitudes of th§ perturbation pressures P, and P_ are given by 

I I c 
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(4.22) 



with 


= (< 5 ^q *^| *S^2 * • • • * • ■ • ) 


$ - ($q . i ^2 *^2 ’ ' " *m * * * * ^ 






■ I **/'’ 


) tan dp-j 


« 3 - k,«, - « 2 f 0 
^0 ’ *k 1 i* r l' /A k 1 » 


cot 2 . sin u, sin 2w, ... sin mu 


7T= (< 

^3 = ^ k 2^1 “ %^ S * 
If = f* U y 2 

Ju h Tk 


,...) 


1 - 1 


2 y 2 


Po tt- 


R 12 h 2 fi.k 2 i 


(p-j)e p 2 dp 2 


ak 2 a. 


< 


■ f 12 *! (p,: 

^R 12 h 2 fi,k 2 s. 


2 P 2 a+ 2 


2 - i ?2 tan a 0 dp. 


‘2 p 2 


1 


*so < r i> r j e ~ 2 2 »: 1/2 

ak 2 l 1 ‘ siku 


(4.23) 


Introducing the expressions for P-j and from equation (4.22) Into 
equation (4.21) we can write the perturbation velocity components of u-j and 
U 2 ^ concise form 


rl 

■ K r/ 

U „1 - K ol^ 

U t 1 ■ K ,l^ 

r2 

• K r2^ 

U <j2 ' K «2^ 

u i2 ‘ K t2^ 


(4.24) 
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where 


K rl = + 


-l^lV 

Z 1 e 

4iryBM 2 


ff a ( z ]| W^x 

U/M f-v 01 h ’*1-*. 


ol 


\1 


— 1 co- Z. /X 1 
1 Z, e 11 1 

♦ — 1 2 l 

4iryBM Z k 


V cot 2 e M 




df-, dz 1 


-Z, 


2r 


Vo I I l 

32iryB M ~ k ? 8. 




exp 1[f A Z, + k,(e - <p - <p )] dc 0 dz 


4 2 2 


2 2 


o2 


v III If* 


32iryB H j k ? i 


3 ^50^3 6Xp ^ f 4 Z 2 + V 0 " v 2 v 2^ df 2 dz 2 


ft tan e. - cot e. 

h„ h 9 

K = — — K ? 

T 2ft - tan e. tan e. 

h l h 2 


(4.25) 


are row vectors. In the above, the Integration over z^ Is an Indefinite 
Integral valid over the range (-«>, + ®), the Integral over Is a definite 
Integral extending over the domain (C 1LE » for the front P r °P el1er and 

over the limits (f 2 LE’ ^2TE^ ^ or rear P r °P e ^ er - 

BOUNDARY CONDITIONS AND BIRNBAUM COEFFICIENTS 
The equations to the upper and lower surfaces of the front and rear pro- 
peller blades can be expressed In terms of the mean camber line and thickness 
and may be written 

V - z il ( V y i> Z L1 - Z L1<W 


Z u2 = Z u2 (r l ,y 2 ,t l ) Z L2 = W’W 


(5.1) 
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assuming that the profile of the blade is independent of the radius r. In 
the coordinate system chosen, the front rotor blades are stationary and the 
equation of the front blade surface is Invariant with time. But, since the 
rear propeller moving relative to the front prop with an angular velocity « r , 
the rear prop surface is a function of time. We can write the dimensionless 
equation for the m-th twisted, tapered and swept blade in the (r, e, z) 
reference frame as front: 

F = z, - r, cot a, tan (e - <p, - 6,1 - - Z' + Y' cot a, + z'c, esc a, = 0 

II 11 lm ul ll 111 1 

F * Z 1 - r l cot “2 ta " ( 9 - *2 - 4 2 m f i^ 1 ) 

- C 02 - Z 2 + Y 2 cot a 2 + Z 2 C 2 CSC a 2~ ° 

(5.2) 

where 

the dimensionless axial displacement of the propeller from the 
origin (fig. 1(b)) . 

Yj, Z j the dimensionless Y and Z coordinated of the blade half-chord 
line S (fig. 1(a)) such that Y' = Y/R and Z' = Z/R. 

C 1 dimensionless blade half-chord at any radius (fig. 2(a)). 

Z^ dimensionless ordinate of the blade profile, 

blade angle distribution. 

tp^ off-set angle of the first blade (fig. 1(a)). 

the subscript 1 = 1,2 to designate the front and rear prop blades respec- 

tively. From these the upper and lower blade surfaces are given by the 
equations : 
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Z 1U = ^01 + l \ - Y i cot a i ~ Z 1U C 1 CSC a l + r i COt c ’ tan ( e - v-, - « 1(n ) 

z 1L = Cq-j * — °] ~ Z ") i_^"| ®i ^ i " 1 1 cot a-j tan (0 - <p-| - 

Z 2U = f 02 + Z 2 ” Y 2 COt a 2 ~ Z 2U C 2 CSC °2 + r l Cot a 2 tan ( e ~ ^2 ' 4 2m + 

Z 2L = f 02 + Z 2 ~ Y 2 cot a 2 ~ Z 2L C 2 CSC a 2 + r l cot a 2 tan ( 6 " <P 2 " 6 2m + 

(5.3) 


The flow tangency condition for the front prop may be expressed as 


r Ul 


!!iiu 

dy( 


All 

J i r 


LI 


= 0 + 


^ik 

dy-j 


'U 


1 z 1 

j iy 


0 - 


(5.4) 


where , and U^ z , are the y.j and z.j components of the resultant 
velocity at the front prop due to both the propellers. For the rear prop, the 
surface boundary condition Is obtained by equating the substantial derivative 

1 

DF/Dt = 0 where F(r, e, z, t) = 0 Is the equation of the rear prop blade 
surface. For an arbitrarily chosen first blade, m = 1, this may be written as 


U 2r s1n 2^e - i 2 ♦ ^ ♦ 


2 / _ 

3U 2o “ 2 < t 2 c 2 + tan a 2 )U 2z C0S l 0 " v 2 + ~ 


nt 


1 


= 2 tan a 2 + 2 (t 2 c 2 + tan c* 2 ) cos I © - <p 2 + - — J - 4« tan © h 


(5.5) 


Since the dimensionless free stream for the rear propeller has the velocity 
components 


U 2r - »rt » U 2e * ta " 6 h, - tan “2 *.V* U 2z ' 1 ‘ “z2,, <*:*> 

The boundary conditions on the rear blade will be satisfied at an arbitrary 
time t-| = 0 on the upper and lower surface and can be written 
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U 2r sin 2(e - (P 2 ) + 2U 2e - 2(t U2 c 2 + tan a 2 )U 2z cos ~ ^ 

2 

= 2 tan a 2 + 2( T u 2 c 2 + * an ^ cos (® “ <P 2 ) _ tan e h 


J 2r sin 2(6 - v 2 ) + 2U 2e " 2 ^ t L 2 C 2 + tan “2^ U 2z cos ( e " <P 2 ) 

2 

= 2 tan <* 2 + 2 ( T |_2 C 2 + tan a 2^ C0S ~ v 2^ " tan 0 h 


(5.7) 


The velocity components U.^, U^, and U^, can be expressed In terms of the 
components (u^, , u^) and (u r2 , u^) as 


U lr ' u rl * u r2 


U, . = cos a, + tan 0 . sin a, + u , cos e, + u , sin e, 
ly‘ 1 h, 1 al 1 1 


+ u o2 cos © 2 + u t2 sin © 2 


U 1z , = sin - tan cos sin cos 

- u sin e„ + u „ cos e„ 

a 2 2 t2 2 

(5.8) 

Substituting Into equation (5.8) we can write the boundary conditions for the 
front prop as 

Ml’ * k ,i *‘5V ♦ ( k « 2 *‘1 > * k , 2 *U > )« 

(RJ 


Val’ * K ,1 5 2 jV ‘ {*Al ] * ^A*)* -*2 ' 


(5.9) 
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where 

am 


are defined by 


and % 


A (1) 

A 11 

= sin e.| + 

T (1) 

T U1 

cos 

e l 

A (1) 

13 

= sin e 2 + 

(1) 

T U1 

cos 

6 2 

A (1) 

«2i 

= sin 6-j ♦ 

(1) 

T L1 

cos 

6 1 

A (1) 

m 23 

= sin e 2 + 

T (1) 

T L1 

cos 

6 2 

r (D 

L 1 

= sin a f - 

- T (1) 

T U1 

cos 

a r 

c (1) 

L 2 

= sin a p - 

- T (1) 

T L1 

cos 

a r 


A (1) 

rt 12 

= - (cos 

6 1 - 

T (1 ) 

T U1 

sin 

A (1) 
A 1 4 

= - (cos 

e 2 ~ 

T (1) 

U1 

sin 

A (1) 

a 22 

= - (cos 

6 1 - 

T (1) 

T L1 

sin 

*0> 

«24 

= - (cos 

e 2 - 

T (1) 

T L1 

sin 

(cos 

a + T (1) 
r T U1 

sin 

a r ) 

tan e 


(cos a r + T<|> sin a f ) tan e 


e l> 


0 2> 


V 


e 2 ) 

h l 


'1 


(5.10) 


Considering the rear propeller, we have on substituting for the velocity 
components (u p2 . u 02 , u z2 ) the following 






U 2r * 

u , + u „ 
rl r2 



1 

26 " U cr 

sin 

e h 

h l 

- “tl 

COS 0. 
h l 

+ u , sin 
o2 

e h - 
h 2 

U^ 2 COS 

e h 

h 2 

22 = U ol 

cos 

e h 

h l 

+ “,1 

sin e hi 

+ u 0 COS 
©2 

e h + 
h 2 

“tZ S,n 

e h 

h 2 


• (5.11) 

i 

Into equation (5.11), the equations for the boundary conditions on the upper 
and lower surface of the rear propeller become 


( K rl S 6Q ) + K 0 1 S 61 ) + K t 1*62*) j * + ( K r2 A 60* + *02*63* 

I 

(k S (1 ^ ♦ K ' A (1) ♦ K S (1) W + U A (1) + K A (1) 
\ rl A 70 ol A 71 t! A 72 r V>2 A 70 o2 A 73 





<s 1 ) 


= <i. 


(S^ 


(5.12) 
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m ( V 

where we have defined A^' and as follows. 


2 - 

A fe1 = 2[s1n 0 h - ^ T U2 C 2 + tan a s^ cos ~ 4> 2 ) C0S 9 h ^ 

2 

A &2 = - 2[cos © h f (ty 2 c 2 + tan a 2 ) cos (e - <p 2 ) sin © h j 


A ?1 = 2 [ s 1 n - ( t l 2 c 2 + tan 01 s ^ cos ” <f> 2 ) cos e h ^ 


2 

A ?2 = - 2[cos + (t l2 c 2 + tan a $ ) cos (e - <p 2 ) sin © h ] 


A 63 = 2[s1n e h - {t M ,c, + tan a t ) cos (0 - tp,) cos e h ] 


'64 


h 2 V ‘U2"2 


2 

2[cos © h +■ (-Ty 2 c 2 + tan c* 2 ) cos (9 - <p 2 ) sin ] 


A ?3 = 2[s1n 6 h - <\2 C 2 + tan a s* cos < 6 ~ ^ C0S 6 h ^ 


A 74 = - 2[cos © h + ( t l 2 c 2 + tan a s > cos (e - tj> 2 ) sin ] 




(SJ 


= 4 tan « 2 - 2(1 + 2Q) tan 0^ + 4 ( T u2 c 2 * tan cos ( e ~ 4> 2 ) 

* 0 ) 


, 6Q - sin 2(6 - <p 2 ) 


(S ) 
<4. 1 


7 = 4 tan a 2 - 2(1 + 2Q) tan * 4 (t |2 c 2 * tan a s ) cos ( 0 - <p 2 ) 


A^ } = sin 2(0 - v 2 ) 


(5.13) 


Therefore, the boundary conditions, equations (5.13) and (5.12) can be written 
In the compact form 
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M + 

L# = 

Y 1 

K 

- 

* K xl S 12 J 

L = 

K o2 A 13^ 

+ K t2 S 14 ) 




m = 

Y 2 

M 

- kAA 

. k t1 a* 2 ' 

N = 

*02*23* 

+ K t2*24* 



p^ + 


y 3 

P = 

kAI } * 

rf * 

K ,1 5 M 

Q = K 

r2 S 60 ) + K o2 

A (1) 

*63 

* \2*M 


S$ = 

y 4 

R = 


* 


S = K 

r2*70* + K o2 

A (1) 

*73 

* K T2 i 13 1 











(5.14) 


The four equations may be written as the pair of matrix equations 

V * L N» * r T 

V* Q^- r 2 

( 5 . 15 ) 

wjhere the matrices K^, L^ t P R , Q s , and r 2 are defined by 



( 5 . 16 ) 

Solving this pair of nonhomogeneous simultaneous equations for A and B we 
obtain the Blrnbaum coefficient vectors A and B as 

^ = ( L N K M “ Q S p r) ( L N r i “ Q s r z) 

® = ( K M L N ~ P R Q s) ( K M r i ' P R O 

( 5 . 17 ) 

The matrices K^, l_ N , P R , and Q<. are primarily aerodynamic In nature whereas 
the matrices r 1 and r 2 are purely geometric representing the blade section 
and blade planform. Truncating the Blrnbaum series at m = M* we have 
(M* + 1) coefficients for each of the front and rear propeller. Considering 
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P* points on each side of the blade at which the boundary conditions are 

satisfied, we have 2P equations to determine the 2(M* + 1) coefficients of 

the front and rear propellers. Hence, we must have 2P* = 2(M* + 1) giving 

the relation between the number of points on the blade and the number of 

coefficients of the Birnbaum series. The matrices K u , L.., P_, and Q_ are 

M N K c> 

of order (2P* by 2P*) and the matrices and r 2 are of order 2P* by 1. 
Each of the vectors and ft is of order 2P* by 1. 

The vectors $i and ft are functions of the off-set angles and ip^. 
For purposes of numerical study, we may set ^ = 0 and <j> 2 = & r t = 

Since the blade passing is periodic, it is only necessary to vary the time t^ 
in the range 0 ^ t^ « and obtain the values of sA and ft at desired 

Intervals. 

From a knowldege of ft and ft the pressure loadings and on 

the front and rear propellers can be calculated using equations (4.4) from 
which the thrust and torque for each blade can be calculated for performance 
estimation of the counter-rotating propeller system. If the free stream were 
steady, «.j = 0. The Inflow into the rear propeller would experience the 
inviscid periodic wake of the front propeller with an angular frequency 

Ca>2 = - 

The above solution is applicable for calculating the unsteady and steady 
aerodynamic behavior of the blades which are assumed to be rigid. However, 
this assumption may be relaxed by allowing the blade geometry parameters 
, Z'y and in equation (5.7) to be functions of the blade loading and 

stiffness characteristics for an aeroelastic analysis of the system. It has 
been mentioned earlier (cf. eq. 3.46) that the time dependent perturbation 
pressure P at any point in the flow equals the sum of the perturbation 
pressures and P^ due to the front and the rear propellers. This 

can be calculated from eq. (3.44). 
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CONCLUSION 

The lifting surface theory for the counter-rotating propeller has been 
presented above. It Is possible to obtain the acoustic modes and the cor- 
responding sound pressure levels of the combination using the Blrnbaum coeffi- 
cients d and <3 together with the Blrnbaum series. Further, the 
m m 

aerodynamic performance of the dual rotation propeller system can also be 
calculated. 
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